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Abstract 

We examine the interplay between the percolation and the deconfinement phase transitions of 
Yang-Mills matter at finite temperature, quark chemical potential /xg and number of colors Nc- We 
find that, whereas the critical Nc for percolation goes down with density, the critical Nc for con- 
finement generally goes up. Because of this, Yang-Mills matter falls into two qualitatively different 
regimes: the "low-A'^c limit", where percolation does not occur because matter deconfines before 
it percolates, and the "high-A'c limit", where there are three distinct phases characterizing Yang- 
Mills matter at finite temperature and density: confined, deconfined and confined but percolating 
matter. The latter can be thought of as the recently conjectured "quarkyonic phase". We attempt 
an estimate of the critical Nc, to see if the percolating phase can occur in our world. We find 
that, while percolation will not occur at normal nuclear density as in the large- A'c limit, a sliver 
of the phase diagram in Nc, energy density and baryonic density where percolation occurs while 
confinement persists is possible. We conclude by speculating on the phenomenological properties 
of such percolating "quarkyonic" matter, and suggest avenues to study it quantitatively and to 
look for it in experiment. 



1 



I. INTRODUCTION 

The "large number of colors" approach has been a promising way to simplify 

some of the tremendous mathematical difficulties inherent in handling non-perturbative 
features of Yang- Mills theory. The idea is to take the number of colors N^. to infinity while 
taking the Yang-Mills coupling constant gyM to zero in such a way that Qym^c = A stays 
constant, defined at some perturbative fixed scale. Numerical results, obtainable by plugging 
in Nc = 3, should be correct within ~ 30% or so, and hence this simplified theory should 
be enough for a qualitative estimate. 

While this theory shares with QCD its non-perturbative nature (strong coupling arises at 
a scale ~ N^, parametrically similar therefore to the QCD scale of Aqcd — 250 MeV), this 
approach has led to some important qualitative results: the fact that in a confined regime 

n n 

mesons are quasi-particles p| while baryons are classical states [2[ can be explained in this 
large-Nc limit. Features of QCD such as the dominance of planar diagrams (and hence the 
string description of gluon propagators and extension into the gauge/string correspondence 
[3I) and the OZI rule are also well explained with A^^c-counting. This has made large a 
useful tool for phenomenological as well as theoretical analysis jj]. 

The large-Nc limit, however, has some qualitative divergences from physical QCD too, 
divergences too big to be put down as a 30% effect. Due to the identification of confinement 
with center symmetry restoration js-?], deconfinement is a first order phase transition in 
the large- iVc limit provided the number of light fiavors Nf <^ N^, it is a smooth crossover in 
our Nc = 3 world |8|,b|. 

Nuclear matter is a tightly bound crystal in the large- A^c limit 
in our world 
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10| , whereas it is a liquid 



ll-ll9|. The latter feature is a consequence of the fact that in the Isnge-Nc 
limit the inter-baryon binding energy scales as the baryon mass, A^^cAqcd- In reality, the 
scale of inter-nuclear forces is around ~ (9(Aqcd/(10 — 100)), a "hierarchy problem" which, 
given the soundness of the large- A^ description, needs to be resolved. 

Given the considerations above, a phase transition in Nc, between A^c = 3 and Ac — ?■ cxd. 



is a plausible resolution of some of these issues |20l-l22[|. The existence of two possibly linked 
transitions in Nc is in fact fairly certain, due to the arguments above: if we could keep 
Nf constant {Nf > 1 for baryons to exist) and increase N^, we would find a critical point 
for confinement at zero quark chemical potential (somewhere between the cross-over in our 



world and a first order transition at Nf/Nc — )■ 23|, |2J|) and a liquid-crystal transition for 
matter at high chemical potential (since the laige-Nc matter is crystalline 10| , and a crystal- 
liquid transition is typically associated to a phase transition due to a change in translational 
symmetries) . 

As discussed in [20I (and hinted at from gauge/string calculations 



251, 



26|). the crystal- 



liquid transition is linked to the nuclear matter hierarchy problem: the classical picture 
of the baryon necessarily entails an Nc much larger than Nn ~ (9(10), where is the 
number of neighbors in a densely packed system. Below this limit, one can not ignore the 
Pauli exclusion principle in the non- color part of the baryonic wavefunction. This raises the 
energy cost of compressing baryonic matter by ~ NcA^^-^, and hence most likely lowers the 
equilibrium density to values lower than ~ ^qcu (i'^ fact, even <^ itl^). Since inter-quark 



interactions are suppressed by the confinement scale, and pionic exchanges are ~ e 



the 



nuclear forces get weaker, hence the critical point of the nuclear liquid-gas phase transition 
happens at T,fiQ < Aqcd- 

Going further is hampered by the fact that our understanding of corrections to Nc in this 
regime is limited. The regime in which this transition occurs is inevitably strongly coupled, 
making perturbative calculations untenable. Standard methods of lattice QCD can not be 
used, since fiq/T ~ O (1) (27H29[|. And, as discussed in [30i], gauge/gravity techniques are 



also unreliable since this transition, by its very nature, is quantum- gravitational, something 
of which we have a very limited understanding. 

The only possible way to move forward, then, is to investigate models which are simple and 
qualitative, yet are universally applicable. One suggested way to describe phase transitions in 
Yang-Mills is via percolation 3l|-l33|. The idea is that, at increasing energies, the increasing 



parton densities will make partons of different hadrons "overlap" as their interaction cross- 
section becomes of the order of inter-parton spacing. It is logical to associate this transition 
to deconfinement, where a quark can propagate throughout the hot medium rather than 
being confined to the hadron size, whose natural scale is ~ ^qcu ~ 1 While these 
analogies might touch on deeper conceptual issues [34|], the percolation picture of confinement 
misses the order of the phase transition both at A^^ = 3 and A^ — ?■ 00, so its direct relevance 
to confinement is questionable. 

There is however a newly conjectured regime where the percolation picture might be 
viable: it is the proposed "quarkyonic matter" at low temperature (below the deconfinement 



N large 

c ° 

i 

III Quark-Gluon plasma 

p~n' 


T ~N° 




c C 


II 


I confined vacuum 


Quarkyonic matter: I 





Nuclear matter: I 
1 = >- 



1/2 



P~N, 



' C / 



y QCD 



N^=3? 



T 

A.. 



baryon 
density 
jump 



y QCD 



FIG. 1. (color online) Larger left panel: The phase diagram in the Nc ^ oo limit, where the 
deconfinement line becomes flat as quark corrections vanish, and a first order phase transition, with 
the baryon density as an order parameter, rises up vertically w.r.t. temperature at fiq ~ Aqcd- The 
scaling of the pressure with Nc is shown in the different phases. Smaller right panel: a qualitative 
sketch of the expected situation at Nc = 3. 



iemperature T^) and moderate density (one baryon per baryonic size, fig = fis/Nc — A 
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qcd; 



37H43| , which is confined (the excitations at the Fermi surface are baryonic) but "quark- 



like", in that pressure and perhaps also entropy density feel the quarks below the Fermi 
surface and consequently scale as N^, as opposed to N^^. 

Unlike in 



31 



33| , the "quarkyonic" transition is thought to be distinct from deconfine- 



44|-|47l| at low temperature but high 



ment, to be found in lower energy heavy ion scans 

baryo-chemical potential = ^cUo £a^ description of deconfinement at finite density as 



percolation was also postulated in 
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49|). 



The motivation is that we expect confinement to disappear when the quark-hole exci- 
tations at the Fermi surface are of the same order as the ant i- screening gluon loop. Since 
quark-hole excitations ~ ^qNcNj while gluon loops ~ N^, and since the strong coupling 
constant becomes non-perturbative at Aqcd, dimensional analysis suggests that confinement 



disappears at a chemical potential fiQ ~ ^JNc/NjAqcj^. 

Thus, the phase diagram at A^^c — ?■ oo looks like the one in Fig. [1] (left) [355: the deconfine- 
ment line becomes infinitely flat. At the same time, the transition to "nuclear matter", with 
the baryonic density as order parameter, becomes infinitely sharp since the baryon mass 



~ A^^cAqcd- Therefore baryons drop out of the confined vacuum partition function entirely, 
but continue to be present at fig > Aqcd- 

Hence, "nuclear matter" at fig ~ Aqcd should, at large N^, be in the confined phase. In 
configuration space, however, inter-quark distance ~ Nc^^'^'- for large enough Nc, then, one 
should be in the confinement regime yet somehow neighboring quarks should be so close 
that asymptotic freedom applies. The authors of [sS^ proposed a solution to this seeming 
contradiction by postulating matter in this regime is "quarkyonic" , with quark-like degrees 
of freedom deep inside the Fermi surface (and hence a scaling ~ A'^^ for the pressure) but 
baryonic excitations on the surface. 

The possibility of discovering the quarkyonic transition assumes, however, that physics at 



high chemical potential is qualitatively the same when is varied from 3 to infinity. In [2l| it 
has been shown that for a wide variety of reasonable propagators at a fixed baryonic number 
density of = Aq^^/S a percolation transition is found as Nc is varied. If one identifies the 
percolation transition with the quarkyonic phase, deconfinement and percolation are indeed 
separate, and they cover different regions not just in T and /ib, but also along Nc. 

In this work, we extend the results of 2l| to variable density and non-zero temperature. 
The purpose of this exercise is to determine the role of percolation in the full T-fi-Nc phase 
diagram, and to see whether percolation is involved in the physical Nc = 3 world, or, instead, 
whether this transition divides our world from the "truly large- A^c" regime. 

Specifically we aim at determining whether there is a region, in the T-fis plane, where a 
percolating yet confined phase is likely at Nc = 3, Aqcd < Pq < a/ Nc/NjAqcu- If so-l] this 
would be the natural region to investigate for quarkyonic effects in experiment. We also ask 
ourselves whether percolation is related to the more usual liquid-gas phase transition, and 
if its onset therefore accounts for the large phenomenological failures of the large- A^c picture 
in this regime 20 1. 

We close with a discussion outlining what an effective theory for percolating matter 
would look like, and suggesting ways of looking for it in both lower energy experiments and 
astrophysical searches (neutron and proto- neutron stars). 

^ In our world of course Nf = 2 if the strange quark is heavy, and Nf = 3 if the strange quark is hght. Since 
the bare strange quark mass ^ Aqcd, it is far from clear which limit applies, yet this is the crucial question 
determining whether Nf /Nc is an expansion parameter at all. The mass of the strange quark might well 
be the crucial qualitative uncertain driv ing factor in our results, and hence fundamentally determining the 



nature of the QCD phase diagram 
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II. THE VARIABLE-DENSITY PERCOLATION MODEL 



The strategy used to investigate the percolati on p roperties of high-density baryonic mat- 
ter is a generahzation of the model presented in 2l|, to which we refer for further introduc- 



tory details. In our description, baryonic matter is arranged in a cubic lattice, with a baryon 
sitting at each lattice site: its quarks will be randomly, independently positioned according 



to a hard-sphere distribution with radius I/Aq^d- /(x) oc 0(1 — Aqcd|x — x'^^"*'^'' | ) . In |21| . 
each sphere touches exactly its six neighbors, that is, the lattice spacing is fixed to 2/Aqcd 
and the density p is therefore fixed (we will relate p to the thermodynamic baryonic density 
ps in Section ITn|) . Replacing the cubic arrangement with another regular 3D lattice would 



have changed the percolation threshold by O (30%) or so 5^, and hence not impacted our 
results qualitatively. 

The generalization to a variable-density setting is realized by the introduction of the 
parameter e, defined as the ratio of the lattice spacing over twice the spheres' radius: thus. 



the density pg examined in 2l| had e = 1. Since at Pq each baryon occupies a volume of 
(2/Aqcd)'^, we now have 

_ 1 A^ 1 
P(e) = Po^ = -f^^ ■ (1) 

In support of this classical, static description of baryonic matter, we keep an eye to the 
large- limit and note that the propagation speed of Fermi-surface baryons. 



1/ a/Ac/A/ in the confined regime, is parametrically smaller than the characteristic mo- 



mentum of quar! 



at larger A^ 



is, ~ A": hence, percolating quarks see the baryons as quasi-static (a crystal 



lOj . and, presumably, a disordered "glass" at smaller Ac). Deviations from 



the "baryons are spheres" assumption might become significant when the number of colors 



approaches the number of neighbors of a densely packed system [20], which in 3D means 
Ac ~ (9(10). This, as we will show, coincides roughly with the percolation threshold for 
"sensible" choices of the parameters at play. 

Let us briefly summarize the results in 2l| before generalizing to the present setting. 



We assume a probability ("squared propagator") for two quarks in different baryons to 
exchange energy/momentum with the essential properties of (a) getting weaker as A/Ac at 
increasing Ac, and (b) dropping quickly to zero around some confinement scale tt/Aqcd oc 
O (1) Ac . These requirements enforce the relevant physics of the problem; we then consider 
two representative choices for the "propagator", inspired by the Gribov-Zwanziger theory 



53l . |54| : the step- function in coordinate space and the step- function in momentum space 
(that is, the squared Fourier transform of a p-space theta-function). These can arise out of 



chromo-field interactions such as in [48 
by 

A y 



49| . Their expression, suitably normahzed, is given 



F.(.)^-e(i- 

FK{y) 



X 2rl 



2 ( y 



(2) 
(3) 



respectively, where y is the inter-quark distance in physical units. The two parameters are 
~ 1 and A ~ 1 ('t Hooft coupling), controlling respectively the range and the intensity 
of the interquark exchanges. 

Conceptually, these definitions leave some ambiguity of what is "propagated" . Our work- 
ing hypothesis is that Eqs. |2] and |3] represent tunneling-driven interactions of quarks from 
different baryons exchanging conserved quantum numbers (spin, flavor, energy-momentum) 
concurrently with some global color-neutralization mechanism. We will leave discussion as 
to how this could happen to Section IIVI We note, however, that similar propagators have 
already been used in the context of quarkyonic matter, having been instrumental in the 
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56| (the results in [37| are based on 



study of the conjectured quarkyonic "chiral spirals" 
a propagator of the form of Eq. [3]) . 

With F{y) and /(x) as input, then, a probability p for the exchange between neighboring 
baryons A and B is computed via 



1 - 



xa - XbI) 



(4) 



In 



2l| we concluded that, in order to meet the expectation p — > 1 for Nc ^ oo, the correct 



choice is a = 2: we interpret it as a cross-baryon "interaction", as opposed to the cross- 
baryon "propagation" associated to a = 1. This constraint, derived in the large-Nc limit, 
would leave the question of what would happen at lower Nc somewhat dubious. The physical 
interpretation given in section [TV] however, makes it clear that a = 2 should not depend on 
Nc'. large-scale color neutralization should keep color localized at super-baryonic scales in 
any confined phase, all that can be propagated being flavor, spin, and energy-momentum. 
Hence, nothing prevents the same quark from interacting an arbitrary number of times with 
neighbouring quarks, resulting in a = 2 for all N^. 



The resulting p is a function of Nc (see 2l|] for details) and can be compared with the 



bond-percolation threshold to determine if large-scale correlations occur or not, identifying 
a critical A*"* where the system starts to percolate. Note that direct exchange between 
non- neighboring baryons is neglected by construction. 

However, in the variable-density setting, when e is small we have substantial overlapping 
between the spheres, and non-negligible contributions from non-nearest-neighbor direct ex- 
change: we need to take into account a set of baryon-to-baryon probabilities Pi, one for 
nearest-neighbors, one for neighbors with relative distance of (1,1,0) spacings, and so on. 
In practice, we considered nine probabilities {pi,i = 0, ■ ■ ■ ,8}, associated to nine "neigh- 
boring classes" (we assume the probability of direct exchange between spheres that are 
further apart can safely be neglected): the corresponding relative distances, in units of the 
interbaryonic distance, read: 

(1, 0, 0), (2, 0, 0), (1, 1,0), (2, 1, 0), (2, 2, 0), (1, 1, 1), (2, 1, 1), (2, 2, 1), (2, 2, 2) . (5) 



Thus, for each choice of the parameters (r-r. A), each "propagator" and each e, a set of 
{pi} can be evaluated numerically with Eq. HI by placing the two centers at the desired 
distance. Representative results are shown in Fig. [2] for the Fx propagator. We examine 
9 values of rr E [0.6:1.4], 6 of A G [0.7:1.2], and 13 of e G [0.8:1.4]: the latter corresponds 
to ranging in density from ~ 0.046 Aq^^^ to 0.244 Aq^^^. We generate data for all integer 
from 2 to 80. The {pi} are calculated numerically by sampling the integrand about 4 ■ 10^ 
times per setup. 

The next step is to translate the set {pi} to a single probability p, which will be compared 
to the 3D cubic-lattice bond-percolation threshold pc — 0.2488. This is accomplished by 



means of Monte Carlo renormalization (see, e.g.. Chapter 4 of 57|): we perform a blocking 
step on a cell of sites, mapping the problem to a super-lattice whose sites are the corners 
of the cell; the corresponding p is evaluated by numerical simulation, with the {pi} as input, 
as the probability that two opposite planes on the cell are connected by a continuous path. 
Since po is the only relevant coupling in the blocking-out flow, we are effectively moving on 
a RG trajectory in the nine- dimensional space of the pi, whose fixed point - the percolation 
critical point - lies on the axis Pi = 0,i > 0. The procedure is illustrated in Fig. |3l In 
practice, we start with an empty lattice and, considering all pairs of sites, we switch 
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FIG. 2. Baryon-to-baryon probabilities for the Fk propagator; the reference case is e = r^^ 
1.00. Left top: pQ for varying e. Right top: po for varying A. Left bottom: po for varying 
bottom: pi at all neighboring-classes considered. 

P 

• f 

FIG. 3. Schematic procedure for the RG-step. The 6 = 3 cell on the left, with all its bond 
probabilities pi, becomes an elementary square of the super-lattice on the right, with the associated 
p as its only, nearest-neighbor, bond probability. 

on the links according to the Pi for the corresponding neighboring class. At the end of 
the process we can have a continuous path connecting the 2; = wall of the cube with 
the z = h — 1 wall: the probability p{pi) for this to happen is computed by repeating the 
operation many times (in our case, 20 thousand times per setup). 

9 
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The exact choice of the "crossing rule" -Ri (that is, we get a 1 if and only if there is a 
connection between two opposite wall in a chosen direction, ignoring what happens along 
the other two directions) is arbitrary: for b large enough, all recipes would lead to the same 
final result. 

In the limit 6 — )■ oo, the crossing probability is a step function, zero before percolation 
and one after percolation. When evaluated exactly at criticality, it assumes a universal 
value depending only on the crossing rule considered. In three dimensions, with the Ri rule. 



p*{b) = Hi + with the limit value Hi = 0.265 58|, and leading correction b^^'^, with 
u ~ 0.8765 the 3D percolation critical exponent. 

In order to check the implied assumption of a 6 "large enough", cell sides from 3 to 7 
lattice steps are employed, and the results compared for stability. This technique is first 
tested on the one-dimensional subspace of nearest-neighbor-only pa — )■ p, where we roughly 
reproduce the scaling to the large-6 limit: yet, the curves Pb{po) have a rather 6-dependent 
value p* at the critical threshold: 

= {0.5825, 0.5226, 0.4788, 0.4488, 0.4293} , for 6 = 3, 4, 5, 6, 7 . (6) 

This leads us to define an effective, 6-dependent critical A'^^, the value at which the full 
renormalization step yields the p* for that cell size (practically found by interpolation). We 
are satisfied with the procedure whenever N* shows a substantial independence of b; indeed, 
this will be the case: for example, at e = 0.8, A = 0.7 and = 0.6, we get for the propagator 
Ft (cf. Fig.HD: 

A/"; = {27.9,27.6,27.9,27.6,27.5} , 6 = 3,4,5,6,7. (7) 

The final result is then obtained by the above population of results. The choice of consider- 
ing bonds between baryons far apart has, however, a drawback: namely, the propagator Fk 
(Eq. |3]), as it is, has local maxima around ?/Aqcd ~ Trrj-, Stttt and so on, which dramatically 
alters the results in a nonphysical way; the coincidence of the secondary peaks with neigh- 
boring positions is of course specific to a cubic lattice, and should not occur in a disordered 
" glass of baryons" . 

More generally, however, the power law "tail" in the propagator makes correlation lengths 
diverge at Nc ~ O (1) for densities arbitrarily close to (9 (Aq^j^)). The implementation 
of such a sharp step propagator in momentum space, however, contradicts our intuition 
of statistical physics: a probe charge interacting with a statistical medium is generally 
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FIG. 4. Determination of the 6-dependent critical N*, for Fj- with e = 0.8, A = 0.7, = 0.6. 
Horizontal lines mark p*. 



exponentially screened by charge- hole pairs, acquiring a screening length of dimension ~ 
(or ~ 1/T for T ~ yu). This is a "universal" feature of interacting systems at equilibrium, 
since any backreaction of a charge on a field, to first order, will give a similar effect (the 
absence of such a screening leads to a divergence in the partition function of the fully 
interacting system, the hydrogenic atom's case is a well-studied example of these classes of 
phenomena |59|, l60|). Since /ig ~ Aqcd in the limit we are using, the propagator should 



be exponentially screened at a similar scale. Then, in the following, we will use a "K" 
propagator altered by an exponential screening: 

2 / H-^QCD 



Fs{y) 



X 24 



sm 



-M\y\ 



(8) 



that is, a times a damping exponential factor whose characteristic length is, in practice, 
fixed to 1/Aqcd- Fs, taking the square root and antitransforming, yields 



g'sik) 



dx\/FJx)e''''' oc i atan 



2A, 



Mr 



/ + 1 



A. 



atan 



2A 



Mn 



QCD / ^_ rr 



A 



1 



QCD 



T \ ^i-QCD 

i.e. a "rounded off" step function in momentum space. 

Following the above methodology, we obtain curves in e, one for each propagator {Ft, 
Fk - later ignored - and Fs) and each (rT',A), that can be translated into curves N*ij>) 
through Eq. [1] 

Examples of such curves are shown in Fig. [51 We note they are well fitted to the form 

A 



Kip) 



ip' + B) 



(9) 
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P'-^^QCD P'A^QOD P'A^QCD 



FIG. 5. Representative results, with some choices of rx,X, for the critical percolation A'^^ as a 
function of the baryon density, coming from using the propagators Ft of Eq. [2] (left), Fx of Eq. [3] 
(center) and Fs of Eq. [8] (right). 



We have shown that the percolation transition hne in p-N^ space is strongly curved. This has 



the potential of dramatically altering the conclusions of [2l|], provided the assumption that 
the relevant density is ~ is relaxed: at a greater density, provided baryons still exist, 

A^* can very well be lowered from O (10) to Nc = 3, making the percolating phase accessible 
to experiment. The existence of baryonic states, of course, implies that the relevant density 
is still in the confined phase. In Section llllt therefore, we examine the interplay between 
the confining and percolating transitions. 



III. THE CONFINEMENT PHASE TRANSITION IN N^-ps-e SPACE 



As pointed out in [35|], if one identifies the density of interest for percolation, pB ~ Ap^^j^/S 
(one baryon per baryon size), with the thermodynamic density, percolation happens firmly 
in the confined phase in the large- limit. The reason for this boils down to comparing the 
quark- hole screening with the gluon-gluon antiscreening at large chemical potential (Fig. |6]): 
confinement is broken when the screening by quark-hole pairs at the Fermi surface (which 
decreases the effective coupling) overpowers ant i- screening by gluon loops (which drives 
the coupling constant above non-perturbative values at momenta ~ Aqcd and ultimately 
causes confinement). As anticipated in the Introduction, this gives an estimate for the low- 



temperature deconfinement point as scaling at least as ~ Aqcd ^/Nc/Nf, with higher loops 
giving a {NjNfY>^/^ scaling |36|. 



We also know that - where the chemical potential for deconfinement is zero - does not 
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FIG. 6. Interplay between anti-screening by gluons (driving confinement) with screening at high 



chemical potential. Since any quark-hole 
corrections will not alter this dependence 



oop ~ NfNc while any gluon loop ~ N^, higher order 



36| 



scale with Nf, Nf.: 

- A,e. ; ,r\T.n - ^ (0) ^eiv? ' - (^) ^civ? (10) 

In this Section, we will use the ideal gas Ansatz for the nuclear liquid to translate these 
estimates into an estimate of how p^"^^ and the energy density e™"^^ depend on N^. We note 
that, even in the case of Van der Waals corrections of order O {Nc), the density ps at a given 
T will change at most with terms of O (1) [20], therefore can be ignored for this calculation. 
We start with the ideal gas formula for the density of a relativistic massive gas of fermionj^ 

Pb = '^^^ {Bi{m,T,fiB) - Bi{m,T,-fiB)) , (11) 
the first term accounting for baryons and the second for antibaryons, with 5l|, |61| 



Si(m,T,/i5)= / r / . ^'"^^ XI = (12) 



exp 



1 



= 2^(-l) m -expi^-^)K2[- 

n=l 

and Kn{x) = e~^'^°^^^ cosh.{nt)dt the modified Bessel function. The degeneracy terms 
Qf and Qs are both somewhat non-trivial: gj counts the total "generalized isospin" states 
accessible allowed for the baryon, and so is 1 for one flavor, Nf{Nf — 1) for a higher number. 

If we ignore excited spin states, the baryon will always have spin 1/2, so that gs{Nc) = 2; 
we will later refine this simplification. 



^ For "baryons" to be fermions, we should limit ourselves to integer, odd Nc- as long as we don't include the 
excited spin states, however, the formulas are valid - and have here been calculated - for generic positive 
values of Nr. 
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From Eq. [10], we parametrize the zero-temperature baryonic chemical potential as 

\T = 0) = N^/'N-'^'A^cu (13) 



,, _ ,,conf/ 



(roughly nI^'^Nj ^^"^ baryons need to overlap), and the baryon mass is m ~ A^^cAqcd: we are 
led to change the variables to 
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-m ; a 



-p ; P - 



1^0 



T 



Nf 



note that at zero-temperature deconfinement we have /3 — )■ oo and 7 
critical for confinement as: 



(14) 



Nf. We write the 



nonf (o ^ A-M _ ^^9f9s .3 R 

Pb [p > AqcdJ - TTTiTr 77372 ^QCD^; 



:(iVc,7,/3) 



(15) 



where 

i^2(iVe,7,/3) 



1 + exp 



/3 v^^^T^ 



y {same with /i — )• — ^u} 



ra=l 



/3 



sinh 



^iV,^/3 i^2(n7/3) 



(16) 



This is the general expression for the baryon density at deconfinement as a function of Nc- 
For very low temperatures (/3 — )■ 00, /i = /iq), the second term in B2 vanishes and the 
exponential at the bottom of the first term is infinity if x > and zero if x < 0, where 



X 



(17) 



We then get easily the low-temperature behavior: 



p^°-f(T = 0) = A 



QCD 



Qtt' 



;9f9s 



3/2 
/ 



67r2 ^^°ArJ/2 



(18) 



At T = 0, ground-state baryons are the only possible hadronic degrees of freedom of the 
system. Hence, one can trivially identify ps with the p of Section [Til and directly compare 
deconfinement with percolation: this is done in Fig. [3 In contrast to what we find in 21|, it 
seems that a confined but percolating density at A^c = 3 is possible. The discussion in Section 
[IT] elucidates what [2l| missed: because of the curvature of the density in ps-Nc space, the 
critical drops very rapidly with density, while the density required for deconfinement 



rises with Nr. However, the conclusion made in 



2l| and |20| still stands in that the densities 
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FIG. 7. Log-log plot of zero-temperature deconfinement curves versus percolation curves in the 
plane of pB — J) and Nc- The deconfinement curves (Eq. [T8]) are calculated for Nf = 1,3 and 
are compared to their large-A'^c asymptotic form. The percolation curves are shown with their 
parametrization from Eq. El 



required for it are well away from normal nuclear density, as in the large- A'c limit. Hence, 
percolation is well distinct from the nuclear liquid-gas phase transition and might not arise 
if the strange quark is "too light" (see the footnote [1]). In order to extend our knowledge to 
nonzero temperature, we choose the simplest parametrization for the deconfinement line in 
the T-fiB plane, i.e. that of a quarter of ellipse whose radii are given by the known points 
T = Te and T = 0: 

. « = (19) 

Considering that the transition line is given by the interplay of the matrix elements shown 
in Fig. |6] this elliptical parametrization is actually physically well-motivated, although it 
misses the ~ T/i interference between screening and anti-screening. 

When raising the temperature, we should also include states of higher spin in our calcu- 
lation: in the large- iY^ limit, a spin flip has a cost of ~ Aqcd/-^c 6^, and there can be up to 
{Nc — l)/2 of themlfl We now assume Nc = 2Q + 1 odd integer, and introduce a sum over spin 
states parametrized hj rj = 0,1, ■■■ ,Q {in this way, neglecting higher-spin states amounts 
to limiting all sums to rj = 0), each carrying its degeneracy [21] + 2): this setup replaces the 
factor Qs of Eq. [TTl Thus we write, for nonzero temperature and including higher spins, the 
^ Admittedly, here we prefer to keep the model simple despite the fact that in our Nc = 3 world the cost of 
flipping, say, a proton into a A+ is around Aqcd and not Aqcd/3. 
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FIG. 8. Top: deconfinement line in the PB-^c plane according to Eg. 1201 Bottom: deconfinement 
line in the e-N^ plane from Eqs. [211 [22l The plots are for Nf = 1 (left) and Nj = 3 (right), and 
for various temperatures spanning <T <Tc (denoted by their 9). 



density as 



(27r)3^3/2 QCD 



/ 



'7=0,1. 




(20) 



exp 



3 Nc 1 

2 - /]v7e 



exp 



3 iVc 1 

2 . /tv; 



This relation is plotted in Fig. |8l top panels. As temperature and A^^^ rise, less and less 
energy density is carried by baryons, since the hadronic degrees of freedom are light mesons 
(of mass < 2Aqcd) carrying no baryonic quantum number, and heavy baryons (~ AcAqcd)- 
The critical energy density e^°"^, neglecting the meson mass, is 



„conf 



,7r 



]\f^ I p' 

■^15 ^ 

16 



conf 



(21) 



with the baryonic contribution given by: 



(22) 



a 



^a^ + Nf + r] 



da 



exp 



+ 



3 Nc 1 



a 



^a^ + Nf + ri 



da 



exp 



3 Nc 1 



In the e-Nc plane, at T^, the deconfinement hne 



^conf 



Ni 



(23) 



^15 " 

(with Tc = 165 MeV or so {sl, 1^), is independent of Nc (a vertical line), with a mixed phase 
in the region 

(24) 



15' 



-r 



At T — 7- (that is, — t- 0), the line is found by solving Eq. [22l and similar to the plane 
shown in Fig. [7]- with all energy carried by baryons - in this case we have e = cb- For 
intermediate temperatures, solving Eq. |22] into Eq. |2T] will give an intermediate solution, 
with a non-trivial approach to the T = case for different N^. see Fig. [HI bottom panels. 
The behavior of the energy density is the reason why all curves "curve" (anticorrelate) for 
low-to- moderate values in N^-ps space: at moderate A*",, (including our A^^ = 3 world) 
baryons carry a non-negligible fraction of the energy density even in the vacuum phase, with 
the flavor and spin degeneracy factors beating thermal suppression. At high Nc, while of 
course baryons continue to carry the baryonic number, mesons carry the bulk of the energy 
density. Since the scaling with Nc of the deconfinement line is very sensitive to where one 
is in e-pB (as per Eq. [19]), this interplay can, at lower Nc, change the Nc-pB correlation of 
confinement into an anticorrelation. 

At T > baryonic density and energy density are not the parameters driving the perco- 
lation phase transition anymore. The way percolation treats baryons in Section [IT] does not 
distinguish between baryons created through chemical potential, whose number is conserved 
on average, and baryons created in pairs, whose number fluctuates. Each baryon can give 
rise to quark tunneling and therefore participate in percolating links. 
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For antibaryons at T > 0, a further subtlety is needed: if, by percolation, we simply 
mean exchange of energy and momentum, then, up to Casimir factors of O (1), antiquarks 
percolate as much as quarks. In this case percolation is driven by the total density of 
baryons and antibaryons, Eq. [2D] with the "— " sign replaced by a "+", which would provide 
the definition of p. If we mean instead the delocalization of the quark wavefunction across 
baryons, as we do in Section HV] then only quarks delocalize. In this case, we should recast 
the deconfinement curve in the 'p-Nc plane, where from Eq. [20] we define: 



P 



2^2 ^3/2 ^^QCD 



; ^QCD ^ j ■ 



(2r/ + 2)aMa 



exp 



3 Nc 1 



ITT 



(25) 



At T = 0, where the antibaryon density is strictly zero, these distinctions are insignificant, 
and the comparison in Fig. [7] suffers no problem. At T > 0, however, in principle we need to 
use the temperature and chemical potential to calculate the relevant non-conserved quantity. 
The final result is shown in Fig. [9] assuming the second case, consistently with the discussion 
m Section [IVJ As can be seen, the deconfinement line on the new axis is quantitatively very 
similar to the upper two panels in Fig. [8] Then, the two cases in the previous paragraph 
yield virtually identical regimes on the phase diagram. 

From Figs. [7] and [9] we see that, in the T-ps-Nc space, there can be three distinct phases: 
confined, deconfined, and confined but percolating. In the latter case, arising at high N^.-, 
the Polyakov loop expectation value is still zero and baryons are still physical degrees of 
freedom but {q{x)q{x')) should not vanish at scales larger than the baryon size due to 
tunneling-driven quark interactions across asymptotically large distances. At low Nc, the 
confinement density is lower than the percolation density. Since the percolation transition 
necessitates baryons as physical states, it therefore does not occur, and quarkyonic phases 
such as in 35|] are not realized. 



The critical allowing for a distinction between deconfinement and percolation is typ- 
ically (9(10°) at 6^ ^ 1 and O (10^) at ~ 1. Hence, the percolation phase at N^. = 3 is 
accessible at T ^ Aqcd, but at nuclear densities about 2-3 times that of the liquid phase. 



which, at least using the scaling of Eg. 



lower in T as that examined in e.g. 38 



Ol are not yet confining. This regime is somewhat 



40 



4l| . making it likely that percolation dynamics 



is relevant in proto-neutron stars (a similar regime to that examined in 



lower energy colliders 



63 



641) rather than 



4J-[47|. In fact, given that at Nj = deconfinement and nuclear 
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FIG. 9. Percolation and deconfinement as a function of p (defined in Eg. I25p . for one flavor (left) 
and three flavors (right). 



matter parametrically coincide, the crucial parameter determining the existence of a con- 
fined but percolating phase, rather than an "early deconfinement" at fig ~ Aqcd, might be 



the strange quark mass [50|, ISl] . 



Given the quantitative roughness of the models considered here, however, these are in 
no way definite conclusions. If anything, these results are much more encouraging for phe- 
nomenology than those of [2]]], which suggested that at Nc = 3 the percolating regime was 
strictly inaccessible. Given the uncertainties illustrated above, we will devote the next two 
sections, IIVI and |Vl to exploring some theoretical and phenomenological aspects of the per- 
colating phase, in view of both giving experimentalists and astrophysicists some insight into 
how this phase could manifest, and sketching what the effective theory of this percolating 
phase might be like. 

We close this section by comparing the percolating phase to the more "usual" nuclear 

19 



matter in the large- A^c limit. One can ask how the phase considered here overlaps with 



the transition between the classical baryonic crystal considered in [10| and normal nuclear 
matter. The equivalence between the percolation transition and the onset of the classical 
regime for baryon dynamics can not be exact, since percolation is insensitive to the number 
of flavors (as long as Nf > 1, required for baryons to exist), while the transition is driven 
not by Nc/Nn ^ 1 but rather Nc/ {NfM^) ^ 1, where is the number of neighbors in a 



closely packed system 



20|. 



However, the interplay between this transition and deconfinement has a similar, albeit 



-1/2 



rather than 



Since the critical Nf. for this 



weaker, dependence on Nf y~ 
transition is of O (10), the two transitions do approximately coincide for Nf ~ O (1). This 
suggests that varying Nc and Nf separately could yield extremely non-trivial dynamics. 
Such an "experiment", of course, is only possible on the lattice, perhaps by applying the 



strong-coupling methods of [19!] to the large- A'^c limit 



23 



24|. 



IV. EFFECTIVE THEORY OF THE PERCOLATING PHASE 

In this work we have used a simple but universal model, motivated by what we know 
about A^c-scaling of the thermodynamics of Yang- Mills theories, to map the interplay of per- 
colation and deconfinement across density, temperature, and number of colors Nc- We found 
a hitherto unexplored percolating phase, where confinement persists (the Polyakov loop ex- 
pectation value vanishes and baryons exist as semiclassical soliton states) but quarks are 
able to propagate to arbitrarily high distances via inter-baryon tunnelin g, w hose probability 



53 



54| . Because of 



is governed by a non-perturbative color-neutralizing propagator a la 0, 
this, perturbative quarks and holes should be able to coexist at momenta ~ Aqcd in the 
background of baryonic "classical" potential wells. 

While, in many ways, our pha se bears similarities to the "quarkyonic phase" conjectured 



m 



35| and explored in 



37 



55 



561], there are also differences: in the quarkyonic phase, 
excitations about the Fermi surface are assumed to be ~ A'^ at all momenta, and hence 
entropy continues to be ~ N^, since sub- Fermi surface states carry no entropy. However the 
Gibbs-Duhem relation, linking energy density e and conserved charge density pb to pressure 
P and entropy density s, 

dP P + e- fiBpB 
'=dT= f ' 
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sub-Fermi ^ j^'^"'^"^^"'^^ sub-Fermi 
quark excited quark quark 

FIG. 10. A typical diagram introducing a temperature dependence on the pressure 

seems to demand a. s ^ Nc scaling in an interacting phase where P ~ Nc- If the pressure 
scales as P ~ Ncf(T), then the only way to avoid entropy density to scale as ~ Nc is to 
have an equation of state strictly of the form 

P = N'j,ifiB,T) + N'j^ifiB) (27) 

{without temperature dependence of sub- Fermi degrees of freedom). At T = this is cer- 
tainly the case, but quark- hole diagrams such as in Fig. [10] will inevitably add ~ NcfifJ^sT) 
terms to the partition function, representing excited quark-hole states. Exciting these states 

— 1/3 

will cost momentum ~ /Nc ^ N^Aqco with a < 1/2, hence is not suppressed in the 
percolating phase. Therefore, perturbations around the Fermi surface could still be col- 
ored, but confinement should be maintained at super-baryonic distances. Provided we come 



up with a physical way to realize such a system, it is a reasonable way o: 



percolation transition demonstrated here with the phase conjectured in 35|. 



identifying the 



A physical analogy arising from condensed matter physics is the metal-insulator tran- 



sition 



ectrons 



the critical point of this transition can be understood as percolation of e 
across the semi-classical potential wells generated by the atoms of the material IgS]. Tun- 
neling probabilities, as well as electron-electron interactions, are what drives this transition. 
The analogy with the picture presented here, with the baryons taking place of the atoms, 
is immediate. In this picture, the low-A^^c confined system can only be an "insulator" , with 
quarks of different hadrons not interacting or propagating. At high Nc, however, a "confined 
conductor" phase is possible, in which the low-energy degrees of freedom are not N^ heavy 
baryons, but ~ Nc quark- hole pairs: the quarks and holes are confined in hadrons but, due 
to tunneling, can not be univocally assigned to a given hadron. The free energy Eigenstates 
Note that this analogy is not perfect, at least because the critical exponents of metal-insulator and 
percolation are different; in percolation only the two-point correlation is relevant, while typically in metal- 
insulator transitions higher order correlation functions play a part. The source terms in Eq. 1321 could be 
used to represent such terms. 
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FIG. 11. An illustration of the "free quark" wavefunctions of the percolating phase. Baryons 
play the role of semiclassical potentials, analogously to atoms in a conductor. Quark wavefunctions 
are delocalized by tunnelling, and their Eigenstates are in superposition, analogously to electrons. 
Below percolation, tunnelling probability diverges for an infinitely large system, so quarks remain 
localized with an O (1) baryonic cluster. 



are therefore superpositions of quasiparticle quark wavefunctions across the whole system, 
with characteristic momentum k ~ Aqcd, in the same way as the free energy Eigenstates 
of electrons in a metal are delocalized: the "free particle" quark, hole (g, h) wavefunctions 
(combining color, flavor and spin) are not a continuum in k but obey the Bloch constraint 



<h(a; + PB'^') = </xWexp 



this forces the spectral function Pq{k) to be of the form 

oo 



(28) 



(29) 



n=0 



where Pg{k) is a Lorentzian-type function. Diagrammatically, this is shown in Fig. [TTl 
where the solid lines represent the semiclassical "mean field" baryon potentials and the 
dashed lines the delocalized quark wavefunctions Even if is "high enough" for pQCD 
quark- hole dynamics to be relevant, this density of states is radically different from that of 
a free thermal quark-gluon plasma at high chemical potential, where the spectral function 
is approximately constant, Pq{k) ~ 

The combination between the asymptotically free nearly massless quarks with a spectral 



function inhomogeneous in momentum space such as Eq. [29] is what ultimately ena 
chiral inhomogeneities found in 



55 



56l |. and also in models such as 



66 



67| and 



48 



es the 



49]. If 



such "conductive" quarks are in the asymptotic freedom regime in some limit, their dynamics 
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FIG. 12. The difference between free-theory and quarkyonic-theory Feynman diagram expansions. 
The form factor F{k) is approximately the Fourier transform of the potential in Fig. [TTl 



can be computed perturbatively by adding form factors to quark propagators. For scattering 
processes (such as quark-hole scattering of Fig. [12] and the virtual excited quark of Fig. [TU]) 
the quark and hole propagators will acquire form factors F{k) (Fig. [T2|) 



h.n-2 



(30) 



n = 2 would describe a 3D quasi-perturbative regime, while n = 4 would be close to the 
Gribov limit described in section |Tll The propagators of the non-abelian degrees of freedom 
in Eq. [32] would be similarly modified. The form factor F{k) would be the Fourier transform 
of the lattice of nucleon mean fields shown in Fig. [TT] . If quark wavefunctions are of the form 
of Eq. [28], then lower Fourier components "^^f^^lx) can be color singlets (note that the lowest 
mode is ~ p]/^), while the higher modes are allowed to be colored, provided a compensating 
mechanism (such as in 68N70|) exists, neutralizing the color perturbations over scales larger 



than A^;',„. Such color compensation must actually exist in the QGP as well to avoid 



■■QCD- 



paradoxes described, for example, in 7l(|; however, such sub-AgcD correlations in the high- 
temperature regime would be negligible for any thermodynamic property of the deconfined 
phase, since the microscopic scale of this system is above confinement ~ 1/ {N^NfT) <^ ^qcu- 
This is not so obvious in a confined but "quarkyonic" phase; we therefore must invent a way 
for P,s^Nc to hold at T < and yet color neutrality be maintained at scales < Aqcd- 
Condensed matter physics gives us another example of how this could work, namely 



spin-charge separation in ID systems l73l | (an effect that does indeed seem to be found 
in Non-abelian gauge theories [7^, ItsI i ): i n a ID interacting fermion chain, spin and charge 
generally separate. If, as suggested in [s^,^, quarkyonic matter is governed by dimensional 
reduction, such separation could provide the neutralizing force: confinement would localize 
the color part of the wavefunction only, while allowing spin and charge to propagate as 
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copies of a color- singlet field. 

While the quantitative development of such a theory is left for future work, we shall 
present a sketch of how this works. We use the approach of 55|]; the effective Lagrangian 
along an infinite percolating chain of quarks will reduce, for A^^^ colors and Nf flavors, to 



•S* — 5'2Ar^. [/icolor] + Sn^ [^flavor] 



(31) 



where /iflavor,coior are separate "flavor" and "color" fields; note that each comes with "redun- 
dant" copies of the other. 



Neglecting co 
inspired form |37 



or neutralization, both S2Nf and Siy^ should have the following WZNW- 



55 



76 



Sk[i] = ktr 



1 



2Atx 



(fx e^"'\r^d, 



-1 



+ Ji 



branch- 



(32) 



The charges coming in and out of each effectively ID chain through the branching of the 
percolation clusters can be represented by source terms Jbranch- At the percolation transition, 
they can be defined by the requirement of the 3D conformal invariance of the system. 

If / (a;, £, d^t) = 1, this Lagrangian reduces to 76N84| . and this will be approximately the 
case for the flavor part. For the color part, however, f {■■■) = g{x)h{i, d^i) can be used as 
a mean field: 

baryons: can be represented by "mean field wells", having the form 

(<7)(x) = Go(l-5^^n(x)e^"^^'') , (33) 
where Gq ~ Nc- This forces any wavefunction for i to be centered around Eq. 



neutralization: can be enforced by making /(x,£, (9^£) = 1 trigger a large background 



field 



68| for color non-singlet states of momentum p > Aqcq. This effective mass 



could behave in a similar way as the color chemical potentials discussed in 



69 



7Q|. 



It can be seen that, with this Ansatz, Nc "flavor" excitations of arbitrary frequency and 
"color" excitations of momentum k ~ Aqcq survive. For large Nc, these will dominate the 
entropy. 



We thus recovered the premises of 35] , with equations such as |2] as a phenomenological 
form of the mean-field compensator f{x,p). We therefore arrive, from the percolating side, 
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at a consistent physical justification of tlie Ansatze used in tlie first part of tlie paper. Wfien 
Nc is below the percolation threshold, Eq. |28] becomes unphysical, because the probability 
of the wavefunction to tunnel more than one baryonic distance is vanishing. Above the 
percolation threshold, quark wavefunctions are assigned not to a baryon, but to all baryons 
in an infinite percolating chain following Eq. Quark-hole excitations will obey an effective 
action given by [311 In thermal equilibrium, the flavor part of the wavefunction should yield 
entropy and pressure ~ Nc even if color is neutral at super-baryonic scales. 



Why has nothing similar been observed in gauge/string duality, and how can this phe- 
nomenon be characterized in such a picture? Since the percolation transition itself is driven 
by Nc, to model the percolation point one would have to include leading-order Qs corrections, 
where gg is the string coupling constant [2l|. The percolating phase, however, is in the Iovj-Qs 
limit, and therefore could in principle be seen by constructions which include baryons in the 
semiclassical gravity limit, such as [ssj. Yet nothing in these works suggests that the 
dense phase is anything different from a "dense nuclear gas". Quark wavefunctions might 
be delocalized, yet this results in no additional degrees of freedom at the level of the entropy 
density and pressure. 



The problem is that the argument in 35] assumes asymptotic freedom. Even imple- 

86l . ISTj l do not have asymptotic freedom but rather 



mentations such as Klebanov-Strassler 
asymptotic M = ^ SYM with large A for "hard" momentum exchange. If this transition will 
appear in gauge/string duality, it will be subleading in a'. The percolating regime occurs 
at low Qs (high iVjO but higher string tension a' (lower A), while the "nuclear matter" phase 



discussed in 



25 



26| happens in the weak limit of both Qs and a' . 



Typically, in gauge/string constructions |25|, [26 



baryons are represented by stacks of 



-D7-branes, with the "nuclear matter" phase being represented as a deformation of a string 



hanging from charged D-branes due to the charge on the brane [25|]. The extra entropy 
scaling of quarkyonic matter must therefore be driven, in the gauge/gravity picture, by the 
appearance of KK modes in such a hanging string. We conjecture, therefore, that quarkyonic 
matter of the type we discuss arises at low and moderate a' . Subleading corrections in 
Qs will give rise to the Nc percolation transition between percolating and non-percolating 



matter 



2l|, and perhaps the baryon quantum-to-classical transition discussed in [20|. 
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V. PHENOMENOLOGY OF THE PERCOLATING PHASE 



Our calculations show that seeing this phase transition in future experiments in our 



3 world 



44j-|47j might be possible, provided low-T, high-p^ regions are accessible. 



This makes it desirable to extend the above discussion and develop some phenomenology for 
' ;he percolating phase. As in Section HV] the quantitative aspect of this is left for future work 
but we can let the analogy with the metal-insulator transition guide us: a universal 



characteristic signature of such a system is the appearance of band gaps in the spectral 
function of charge carriers due to Eq. [29] and the (weak) interactions between neighboring 
charge carriers 89 1. 

Such band gaps, of mass ~ 250-400 MeV (well below any resonance mass) could be di- 
rectly detected in electromagnetic probes (the spectral function of e~^e~ pairs in heavy ion 
collisions): as the sketch in Fig. [T3] shows, common q hole — ?• 7 — > l^l^ scattering will give 
an approximately flat spectral function for an unperturbed high-p QGP. If quark wavefunc- 



tions are delocalized across potential wells of size ~ Aqcd, M'^ ~ p^^^ will be suppressed, 
analogously to the scattering of x-rays by electrons in a conducting metal, due to the sup- 
pression of delocalized quark states around that frequency: the form factors of Eq. [30] will 

- — ' 1/3 

be F(k = ) ^ 1, and that will depress the scattering cross-section shown in Fig. [13] 

Below that frequency, color-neutralizing effects might suppress the color-part of the spectral 
function, but the flavor part of the spectral function can still contribute. 

Heavy-ion collisions at SPS and RHIC energies have yielded a continuum reminiscent of 



90 



9l| on the top of peaks associated to the decay of hadronic 



the QGP spectral function 
resonances {u, p,ri, . . .), so perhaps the band gap structure can be searched for in upcoming 



lower energy experiments 



44 



47|. 



Alternatively, quarkyonic percolation as described here could be detected in the phe- 



nomenology of neutron stars and proto-neutron stars [63], l6J]. Quarkyonic matter would 
appear at a pressure about ~ A^^c = 3 times that of nuclear matter at the same density, tem- 
perature and chemical potential, while maintaining a heat capacity and an energy density 
comparable to that of nuclear matter. The extra boost in pressure is analogous to the way 
the electron gas dominates pressure in a metal. Such stiffer equations of state are desirable 

33 



for stars such as 



We would like to thank Irina Sagert for discussions regarding this topic. 
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FIG. 13. (color online) Sketch of what the i'^i spectral function could look like in percolating 

quarkyonic matter, in a QGP and in a resonance-dominated hadron gas. The gap in quarkyonic 

1 /3 

matter arises because no delocalized free quark states can exist around k ^ (the form factor in 
Eg. 1301 vanishes). At frequencies below the gap the color part of the wavefunction stops contributing, 
but the flavor part might still be present. 
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Furthermore, a quarkyonic phase in j)roto-neutron stars might be crucial in the dynam- 
ics of supernovae. The effective stiffening of the equation of state might affect the early 
postbounce supernova dynamics and/or black hole formation times during the core-collapse 
of massive stars. The first effect is interesting in connection to the shock- stalling problem 



found in e.g. 93|. In 63|, |6J] this problem was solved with a more traditional deconfinement 
transition making the equation of state softer, since both equilibrium energy density and 
pressure increase at deconfinement, and the mixed phase drives the speed of sound to zero. 

The percolation transition looks significantly different in a way that might make it easier 
to maintain a shock-wave. In the ^ oo limit the phase transition line is vertical in 
the T-^B plane (regions I and II in Fig. [H The real Nc = "i world has a ~ 30% curvature 
correction), and hence both sides in the Clausius-Clayperon equation diverge. Percolation, 
however, implies a second order phase transition, hence the change in pressure {P/{Tpb) 
jumps by ~ A'^c = 3 when ps crosses the percolation threshold, as quarks start exerting 
pressure) can only be gradual with density (the "jump" is a rapid but smooth cross-over 
at any finite Nc), and there is no mixed phase, or jump in energy density or heat capacity 
(e/ {Tps) stays approximately constant as is varied). It would be very interesting to assess 
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the effect of an equation of state with such a transition in calculations such as |63l.l64j.l94|-l96|. 

In conclusion, we have studied the interplay between percolation and deconfinement in 
Yang-Mills matter at finite number of colors, temperature and density. We find that these 
transitions exhibit a non-trivial dependence on N^-, suggesting that, at least for thermody- 
namics, we can not automatically assume QCD is in the "large- A^^c limit" at T,fiQ ~ Aqcd- 
Our calculations, however, show that the percolating phase could appear for ~ (0.125 
- 3)Aq(^q, provided quarks at this density are still confined. Naive scaling in number of 
colors and flavors suggests they are, although we can not say this with certainty. We have 
speculated what the dynamics of the percolating phase looks like and how it is related to 
popular approaches (such as the gauge/string duality) for describing Yang- Mills matter in 
the same regime. Even if the findings here will not be confirmed experimentally, characteriz- 
ing them more rigorously on the lattice and in the gauge/string correspondence opens quite 
a few questions, the solutions of which could help us clarifying the qualitative structure of 
Yang-Mills theories. As a contact with our A",, = 3 world can not be excluded, we have 
closed by suggesting experimental and astrophysical signatures of the percolating phase. 
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